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ABSTRACT 

Network applications, such as multimedia streaming and 
video conferencing, impose growing requirements over Qual¬ 
ity of Service (QoS), including bandwidth, delay, jitter, etc. 
Meanwhile, networks are expected to be load-balanced, energy- 
efficient, and resilient to some degree of failures. It is ob¬ 
served that the above requirements could be better met 
with multiple disjoint QoS paths than a single one. Let 
G = {V, E) be a digraph with nonnegative integral cost 
and delay on every edge, s, t £ V he two specified vertices, 
and D € Zj be a delay bound (or some other constraint), 
the k Disjoint Restricted Shortest Path (kRSP) Problem is 
computing k disjoint paths between s and t with total cost 
minimized and total delay bounded by D. Few efficient al¬ 
gorithms have been developed because of the hardness of the 
problem. 

In this paper, we propose efficient algorithms with prov¬ 
able performance guarantees for the fcRSP problem. We first 
present a pseudo-polynomial-time approximation algorithm 
with a bifactor approximation ratio of ( 1 , 2 ), then improve 
the algorithm to polynomial time with a bifactor ratio of 
(1 -I- e, 2 -I- e) for any fixed e > 0, which is better than the 
current best approximation ratio (0(1 -I- 7 ), 0(1 -I- -)}) for 
any fixed 7 > 0 [18]. To the best of our knowledge, this is 
the hrst constant-factor algorithm that almost strictly obeys 
the constraint for the fcRSP problem. 
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algorithm, auxiliary graph, cycle cancellation. 

1. INTRODUCTION 
I.l Background 

The disjoint quality of service (QoS) path problem is a 
generalization of the shortest QoS path problem and has 
broad applications in networking, data transmission, etc. In 
data networks, many applications, such as video streaming, 
video conferencing, and on-demand video delivery, have sev¬ 
eral QoS requirements, which require the routing between 
source and destination nodes to simultaneously satisfy sev¬ 
eral QoS constraints, such as bandwidth and delay. In those 
applications, a single link might not provide adequate band¬ 
width, and multiple disjoint QoS paths are often necessary. 
Given cost and delay as QoS constraints, the fc disjoint QoS 
path problem can be dehned as follows. 

Definition 1. ('The fc disjoint QoS path problem) Given 
a digraph G — (V) E), a pair of distinct vertices s, t £ V, 
a cost function c : E ^ Zq , a delay function d : E ^ Zq , 
and a given delay bound D £ Zq , the fc disjoint QoS path 
problem is to compute fc disjoint st-paths Pi,...,Pk, i.e., 
E{Pi) n E{Pj) = 0 for every i ^ j £ {1, ..., fc}, such that 
d{Pi) < D for each i = 1 ,..., fc, and the total cost of the fc 
disjoint paths is minimized. 

This problem is VP-hard even when all edges of G has a 
cost of zero [16]. The hardness result indicates that it is 
impossible to develop exact or polynomial-time approxima¬ 
tion algorithms that strictly obey the delay constraint for 
the fc disjoint QoS paths problem unless P=NP. An alterna¬ 
tive method is to compute fc disjoint paths with total cost 
minimized and a total delay bounded by D (equal to kD in 
Dehnition 1), and then route the packages via the fc paths 
according to their urgency priority, i.e., routing urgent pack¬ 
ages via paths of low delay whilst deferrable ones via paths 
of high delay. Then, the fc disjoint Restricted Shortest Path 
(fcRSP) problem arises as in the following: 

Definition 2. ('The fc disjoint Restricted Shortest Path 



problem, kRSP) Given a digraph G = {V, E), a pair of dis¬ 
tinct vertices s, t £ V, a cost function c \ E ^ Zq", a delay 
function d ■. E ^ Zq", and a delay bound D € Zj, the k 
(edge) disjoint Restricted Shortest Paths (kRSP) problem is 
to calculate k disjoint st-paths P\, ..., Pk, such that E{Pi)n 
E(Pj) = 0 for any i ^ j £ {1, ..., k}, Jf, d{Pi) < D, 

and the cost of the k paths is minimized. 

Previous works on the fcRSP problem are mainly bifactor 
approximation algorithms. An algorithm ^ is a bifactor 
{a, /3)-approximation algorithm for the fcRSP problem if and 
only if for every instance of fcRSP, A runs in polynomial time 
and outputs k disjoint st-paths, and the total delay and the 
total cost of the computed k disjoint paths are bounded by 
aD and (30opt, respectively, where Copt is the cost of an 
optimal solution to fcRSP, a and /3 are positive constants. 
Note that a single factor /3-approximation is identical to a 
bifactor (1, /3)-approximation, and we will use them inter¬ 
changeably in the rest of the paper. 

The fcRSP problem is a fundamental problem in multipath 
routing that has a variety of benefits including fault toler¬ 
ance, increased bandwidth, load balance and etc. Although 
multipath routing is not widely deployed in practice because 
of the difficulty of developing efficient algorithms, research 
interest focus on multipath routing is growing. The reason 
is the development of software defined networking (SDN), 
which seems a wonderful place to deploy multipath routing. 
Controllers therein have global information of the network 
and stronger computational ability, which make it possible 
to deploy some complicated routing algorithm into the net¬ 
work. Our results show that fcRSP admits efficient approxi¬ 
mation algorithm theoretically. This may help to boost the 
development of multipath routing in SDN. 

fcRSP has theoretical values beyond multipath routing. 
It is an interesting bicriteria optimization problem, because 
the shortest path problem is a fundamental problem in the 
area of combinatorial optimization. In addition, our method 
on fcRSP might be applied to other related problems. 

1.2 Related Work 

In general, fcRSP is a budgeted optimization problem: 
given the maximum delay constraint, find fc disjoint paths 
of minimum cost, which has been well-investigated recently. 
Efficient algorithms have been developed for budgeted match¬ 
ing and budgeted matroid intersection [3E]. However, those 
methods can not be adopted to solve fcRSP because fcRSP 
cannot be modeled as matching or matroid intersection. An¬ 
other budgeted network design problem, the shallow-light 
Steiner tree (SLST) problem, is attracting lots of research 
interest. SLST is to compute a minimum cost tree spanning 
a set of given terminals, such that the cost of the computed 
tree is minimized and the delay from a specified terminal 
to every other terminal is not larger than D. This problem 
can not be approximated better than factor ( 1 , ylog^ n) for 
some fixed 7 > 0 unless NP C DTIME{ri°^^°^^) [19j . 
Further, although the inapproximability result doesn’t ex¬ 
clude the existence of polylogarithmic factor approximation 
algorithms for SLST, to the best of our knowledge, no such 
algorithms within polynomial time complexity have been 
developed. The algorithm with the best ratio is a long 
standing result by Charikar et al, which is a polylogarith¬ 
mic approximation algorithm that runs in quasi-polynomial 
time, i.e., a factor- 0 (log^ t) approximation algorithm within 


time complexity Due to the difficulty in de¬ 

signing single factor approximation algorithms, bifactor ap¬ 
proximation algorithms have been investigated. Hajiaghayi 
et al presented an (0(log^ t), 0{log^ t))-approximation algo¬ 
rithm that runs in polynomial time m- Besides, Kapoor 
and Sarwat designed an approximation algorithm with bi¬ 
factor 0(jff^))i where p is an input parameter 

|14 |. Their algorithm is an approximation algorithm that im¬ 
proves the cost of the tree, and is with bifactor (Oft), 0(1)) 
when p = t [14]. It is even more interesting that, for the spe¬ 
cial case S = V, Charikar et al’s ratio of 0(log^ n) [5] (with 
quasi-polynomial time complexity is still the best 

single factor ratio. In addition, the existing bifactor ap¬ 
proximations for SLST are also the known currently best 
approximations for the special case. 

Special cases of fcRSP have also been studied. When 
the delay constraint is removed, this problem is reduced to 
the min-sum disjoint path problem of calculating fc disjoint 
paths with the total cost minimized. This problem is known 
polynomially solvable m- Moreover, when fc = 1, the prob¬ 
lem reduces to the single restricted shortest path (RSP) 
problem, which is known as a basic QoS routing problem 
[3 and admits fully polynomial-time approximation scheme 
(FPTAS) (3(17]. The (single) QoS path problem with mul¬ 
tiple constraints is still attracting considerable research in¬ 
terests. Xue et al. recently proposed a (l-l-e)-approximation 
algorithm [23]. When the delay constraint is on each single 
path, and the cost on every edge is zero, the fcRSP prob¬ 
lem reduces to the length-bounded disjoint path problem of 
finding two disjoint paths with the length of each path con¬ 
strained by a given bound. This problem is a variant of the 
Min-Max problem, which is to find two disjoint paths with 
the length of the longer one minimized. Both problems are 
known to be AP-complete m with the best possible approx¬ 
imation ratio of 2 in digraphs m, which can be achieved 
by applying the algorithm for the min-sum problem in |20l 
121] . In contrast, the min-min problem of finding two dis¬ 
joint paths with the length of the shorter path minimized 
is also NP-complete but doesn’t admit fc-approximation for 
any fc > 1 [3] 1111122] . The problem remains NP-complete 
and admits no polynomial-time approximation scheme in 
planar digraphs m- 

A closely related problem, the fc disjoint bi-constrained 
path problem (fcBCP), targets fc disjoint st-paths that sat¬ 
isfy the given cost constraint c(Pi) < C and the given 

1=1,...,A: 

delay constraint J]) d(Pi) < D. For the fcBCP problem, 

1=1,...,fc 

approximation algorithms with a bifactor approximation ra¬ 
tio of (1 -(- /3, max{2, 1 -|- In^) or a single factor ratio of 
O(lnn) (i.e. bifactor approximation ratio (1, O(lnn))) have 
been developed for general fc in [12], where /3 > 0 is any 
fixed positive real number. Apparently, fcBCP is a weaker 
version of fcRSP, and hence all approximations of fcRSP can 
be adopted to solve fcBCP, but not the other way around. 

The fcRSP problem itself has attracted considerable re¬ 
search interest of a number of computer scientists. [4] and 
[TS] have achieved bifactor ratios of (H--^, r(l+^il2£I±li)(i_(_ 
e)) and (1 -|- y , 1 -I- r) for k — 2, respectively. Based on LP- 
rounding technology, an approximation with bifactor ratio 
(2, 2) has been developed in [9]. To the best of our knowl¬ 
edge, however, no algorithm has achieved a constant single 
factor approximation ratio for decades. Our algorithm is the 




first one showing that fcRSP admits a almost constant single 
factor approximation ratio theoretically. 

1.3 Our Results 

The main results of this paper are summarized in this 
section with proofs detailed in latter sections. 

Lemma 3. The kRSP problem admits an approximation 
algorithm with an approximation ratio o/(l, 2 ) and runtime 
0{{^c{e))^^d{e){n^'^Co%TDL)), where C is the cost of 
an optimal solution and L is the maximum length of the 
input. 

We note that the algorithm is with pseudo-polynomial time 
complexity, because the above formula of the time complex¬ 
ity involves Copt, D, and edge costs and edge delays of the 
graph. However, by applying the traditional technique for 
polynomial time approximation scheme design as in [7], we 
can immediately obtain a polynomial time algorithm with 
bifactor approximation ratio (1 -I- ei, 2 -|- £ 2 ). Below is the 
main idea of the technique; For any constant ei, £2 > 0, set 

the delay and cost of every edge e to 

respectively. Then the time complexity of the algorithm be¬ 
comes polynomial time, since 

0((E E d{e){n^-^Ct,%TDL)) 

<o(t 
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Theorem 4. For any constants £ 1 , £2 > 0, the kRSP 
problem admits a polynomial time algorithm with an approx¬ 
imation ratio (1 -I- £ 1 , 2 -|- £ 2 ). 

The time complexity of the algorithm can be analyzed bet¬ 
ter. Anyhow, to the best of our knowledge, this is the first 
constant factor approximation algorithm for the fcRSP prob¬ 
lem that runs in polynomial time, and almost strictly obeys 
the delay constraint. Note that we can set £1 = £2 = £, and 
the approximation ratio will be (1 -f £, 2-\- e), as we claimed 
in the abstract. 

Our pseudo polynomial time approximation algorithm of 
Lemma[3]mainly consists of two phases. The first phase is to 
compute a not-too-bad solution for the fcRSP problem by a 
simple LP-rounding algorithm as in [^, whose performance 
guarantee has been shown as in the following: 

Lemma 5. The kRSP problem admits an algorithm, such 
that for any of its output solutions, there exists a real number 
0 < a < 2 such that the delay-sum and the cost-sum of the 
solution are bounded by aD and {2 — a)*CoPT, respectively. 

Note that a differs for different instances, i.e., the algorithm 
might return a solution with cost 2Copt and delay 0 for 
some instances, whilst a solution with cost 0 and delay 2D 
for other instances. So the bifactor approximation ratio for 
the algorithm is actually (2,2). The second phase of our 
algorithm, as the main task of this paper, is to improve 
the approximation ratio (2, 2) to (1 -I- e, 2 -|- e) in pseudo 
polynomial time. 


The main techniques involved in the second phase include 
the cycle cancellation method, LP-rounding, and some graph 
transformations techniques. We would like to remark here 
that the cycle cancellation method is a classic framework 
that has been used to solve numerous problems related to 
shortest path, minimum cost flow and etc [mmni]. How¬ 
ever, to the best of our knowledge, in previous works the 
cycles used are computed in a graph which allows either 
negative-cost edges or negative-delay edges, but not both. 
Because no polynomial algorithm is known for computing a 
best cycle for cycle cancellation in a graph allowing both neg¬ 
ative cost and delay. This paper not only gives an approxi¬ 
mation algorithm for the fcRSP problem, but also enhances 
the cycle cancellation method by giving a novel algorithm 
of computing bicameral cycle, which is a good-enough cycle 
for cycle cancellation in a graph where both negative-cost 
edges and negative-delay edges are allowed (The formal def¬ 
inition of bicameral cycles will be given later). That is, our 
algorithm can hopefully improve the approximation ratios 
of other bicriteria optimization problems (or namely bud¬ 
geted optimization problems) in which the cycle cancellation 
method can be suitably applied. 

The remainder of the paper is organized as below: Section 
2 introduces cycle cancellation. Section 3 gives the definition 
of bicameral cycles, and then shows that using bicameral 
cycles in cycle cancellation can result in a desired approx¬ 
imation ratio (1 -I- £ 1 , 2 + £ 2 ). This is the first tricky task 
of this paper. Section 4 gives the algorithm that actually 
computes a bicameral cycle in polynomial time, which is 
the most tricky task of this paper. Section 5 concludes this 
paper. 

2. THE CYCLE CANCELLATION METHOD 
FOR aRSP 

This section will give the key idea of an improved approx¬ 
imation for fcRSP based on the cycle cancellation method. 
We shall first state our version of the cycle cancellation 
method (and point out the difference comparing to previ¬ 
ous versions), then give the improved algorithm that is an 
enhancement to the cycle cancellation technique. 

2.1 The Cycle Cancellation Method 

We would like to start from some definitions and nota¬ 
tions. Let El and E 2 be two set of edges. Then Ei © E 2 
denotes the edge set i?i U i ?2 \ {e{u, w)|{e(u, v), e(v, u)} C 
El U E 2 }, i.e., El U E 2 except pairs of parallel edges in 
opposite direction therein. Let P be a st-path and P = 
{e'{v,u)\e{u,v) € P}, i.e., P is P but with the direction 
of each edge reversed. Moreover, the cost and delay of the 
edges of P are negatived, i.e., c{e'{v, u)) = —c(e(u, v)) and 
d{e'{v, u)) = —d{e{u, v)) for every edge e'{v, u) € P. 

Definition 6. (Residual graph) A residual graph G = 
Gres(Pi, ..., Pfc), with respect to G and Pi, ... , P^ <Z G, is 
graph GU{ U E{Pi))\ U E{Pi), i.e., graph G with 

1 = 1, k 2=1, ...,k 

the direction of edges of Pi, ..., Pk reversed, and their cost 
and delay negatived. Q 

For notation briefness, we use G instead of Gres {Pi, ..., Pk) 

^Note that G can contain pairs of parallel edges in the same 
direction with different costs and delays. Thus, G might be 
a multigraph. 


















while no confusion arises. In literature [12] and [18], the 
authors also employed the cycle cancellation method. How¬ 
ever, in their constructed residual graph, the cost of reversed 
edge e'(v,u) is sat c(e'(v,u)) = 0 instead of c(e'(v,u)) = 
—c(e(u, v)), such that the edges in their residual graphs are 
with nonnegative cost. Then, the minimum-mean-cycle al¬ 
gorithm can be applied therein, and hence a best cycle for 
cycle cancellation, i.e., O with minimized, can be com¬ 
puted in polynomial time m- Note that for a more com¬ 
plicated residual graph as in Definition [d] the method above 
can not be applied to compute a best cycle any longer. 

The cycle cancellation method is based on the following 
proposition which can be derived from the flow theory: 

Proposition 7. Let Pi, ..., Pk be k disjoint st-paths in 
G and Oi, ..., Oh be a set of edge-disjoint cycles in residual 
graph G. Then {Pi, ..., P*,} © {Oi, ..., Oh} contains k 
disjoint st-paths in G. 

Proof. Following the flow theory, we need only to show 
that {Pi, ..., Pk} © {Oi, ..., 0/i}contains an integral st- 
flow of value k and the flow is a valid flow in G. 

Firstly, {Pi, ..., Pfe}©{Oi, ..., Oh} contains an integral 
flow of value k that goes from s to t. That is because s and t 
are respectively with degree k and — fc in Pi U ... UP* whilst 
every vertex are with degree 0 in Oi, ..., Oh ■ So s and t are 
respectively with degree k and —k while every other vertex 
is with degree 0 in {Pi, ..., Pk} © {Oi, . .., Oh}- 

Secondly, {Pi, ..., Pj;} © {Oi, ..., O/,} is a valid flow in 
G. On one hand, every edge appears at most once in the 
flow of{Pi, ..., Pfc}©{Oi, ..., Oh} (i.e. satisfies the capac¬ 
ity constraint). That is because every edge appears at most 
once in {Pi, ..., P*,} © {Oi, ..., Oh}, since{Pi, ..., P*,} 
and {Oi, ..., Oh} shares no common edges, and the paths 
of {Pi, ..., Pk} are mutually edge-disjoint and so are the 
cycles of {Oi, ..., Oh}- On the other hand, {Pi, ..., Pji}© 
{Oi, ..., Oh} contains only edges of G, since every edge 
of {Oi, ..., Oh} \ G has a parallel opposite counterpart in 

{Pi,---,Pk}- □ 

Intuitively, {Pi, ..., Pfe} © {Oi, ..., Oii} is to replace the 
edges of Pi, ..., Pk, which have parallel opposite counter¬ 
part _in {Oi, ..., Oh}, by the edges of {Oi, ..., Oh} \ {Pi, 
..., Pi:}. Such an edge replacement with cycles is called 
cycle cancellation. 

2.2 Cycle Cancellation with Bicameral Cycles 
for fcRSP 

This subsection shall apply the cycle cancellation method 
to improve a solution of the fcRSP problem. Let Pi, . .., Pk 
be a solution resulting from the first phase with large delay, 
say X^iLi d{Pi) > D. According to Proposition]?] if we could 
find a set of edge-disjoint negative delay cycles Oi, ..., Oh, 
then we could decrease the delay of the solution to fcRSP by 
using {Pi, ..., Pj;} © {Oi, ..., Oh} as the k disjoint paths. 
Naturally, a question arises: 

Does there always exist cycles with negative delay in G 
when X]^=i d{Pi) > D holds? 

To answer this question, the following proposition is nec¬ 
essary: 

Proposition 8 . Let Pf, . -., P^ be a minimum cost so¬ 
lution to the kRSP problem that satisfies the delay constraint, 
and Pi, - - -, Pk be k disjoint st-paths- Then {P?, ..., Pfc }© 
{Pi, .. -, Pk} is exactly a set of edge-disjoint cycles- 


The key idea to prove the proposition is that: Every vertex 
in graph {P*, ..., Pj)} © {Pi, ..., Pfc} is with degree 0, so 
{Pi*, ..., Pfc } © {Pi, ..., Pfc} is composed by exactly a set 
of cycles. We note that a generalized version (for network 
flow) of Proposition [H] has appear in [I], so we omit the 
detailed proof here. Let {Oi,..., Oh} be the set of cycles 
of {Pi*, ..., Pfc*} © {A, ..., A}. Note that Eti d{Pi) > 

Eti diPt), so Eti dm = Eti rf(A‘) -Eti d{Pi) < 0 
holds, and hence at least one cycle of {Oi, - - - ,Oh} is with 
negative delay. Hence, the previous question has a positive 
answer, formally as below: 

Lemma 9. If E?=i d{Pi) > D, then there exists at least 
one cycle with negative delay in G- 

From Proposition [7] and Proposition [ 8 ] an idea for the im¬ 
proved algorithm is to find the set of cycles {Oi, - - - ,Oh} 
and use them to improve the k disjoint paths Pi, ... , Pfc 
to an optimal solution P*, ..., Pj}- However, it is hard to 
identify all the cycles {Oi, - - -, Oh} exactly. So the key idea 
of the improved algorithm is generally a greedy approach, 
which repeats computing a “best” cycle (i.e. cycle O with op¬ 
timum ^^^), and using it to improve the current k disjoint 
paths towards an optimal solution, until a desired solution 
is obtained. However, a “best” cycle is probably iVP-hard to 
compute. So our algorithm use a bicameral cycle instead. 
Let Pi, ..., Pfc be the current solution to the fcRSP prob¬ 
lem. Then the main steps of our algorithm roughly proceeds 
as below: 

While E*Li d.{Pi) > D do 

Compute a bicameral cycle O; 

Set {Pi, ..., Pfc} ^ {Pi, ..., Pfc}©0. 

3. PROOF OF LEMMA 3 

This section will first give the formal definition of bicam¬ 
eral cycles as well as the formal layout of our algorithm, 
then show that, by using the defined bicameral cycles, our 
algorithm can output a solution with a bifactor ratio ( 1 , 2 ), 
and terminate in pseudo polynomial time. 

3.1 The Improved Algorithm 

Let Pi, ..., Pfc be a feasible solution to fcRSP. W.l.o.g., 
assume E?=i d{Pi) > D- Let AP — D — Ei=i d{Pi) and 
AC = Copt — E^=i c(Pi)i where Copt is the cost of an 
optimal solution to fcRSP. Let O be a cycle in the residual 
graph wrt Pi, ..., Pfc. Then because the meaning of 
depends on positive or negative d(0) (or c(0)). So we have 
three types of bicameral cycles accordingly, as in the follow¬ 
ing: 

Definition 10. {bicameral cycles} 

1- If d{0) < 0 and c(0) < 0 or d{0) < 0 and c{0) < 0, 
then O is a (type-0) bicameral cycle; 

2- Otherwise, 

(a) If d{0) < 0 and 0 < c(0) < Copt, then O is 

a (type-1) bicameral cycle iff where 

AD = D- Eti diPi) AC = Eti (^{Pi) - 
Etic(A); 

(b) If d{0) > 0 and —Copt < ciO) < 0, then O is a 
(type-2) bicameral cycle iff 


Algorithm 1 An improved algorithm for fcRSP. 

Input: Graph G, specified vertices s and t, a cost function 
c(e) and a delay function d{e) on edge e, and a delay bound 
D G Zj; 

Output: {Pi, ... ,Pk}, the set of edges that compose k dis¬ 
joint st-paths. 

1. Compute k disjoint st-paths Pi,...,Pk by the LP- 
rounding algorithm as in [9]; 

2. While YlLi d-iPi) > D do 

(a) If there exists no bicameral cycle in G, then re¬ 
turn “infeasible”; /*The instance is infeasible.*/ 

(b) Compute a bicameral cycle O; 

(c) {Pi,...,PJ ^ {Pi,...,PJ©0; 

3. Return {Pi,..., Pfe}. 


Clearly, a type-0 bicameral cycle is the first choice for our 
algorithm, because it can be used to decrease the delay of 
the k disjoint paths without any cost increment according 
to our algorithm. Differently, using a type-1 (2) cycle O 
for cycle cancellation will decrease delay (cost), but increase 
the cost (delay) of the current solution. Since computing a 
bicameral cycle is one of two tricky tasks of this paper, the 
algorithm will be deferred to Section 4. The formal layout 
of our algorithm is as in Algorithm [T] 

We would like to give some remarks on the definition of 
bicameral cycles before presenting the ratio proof. Firstly, 
a bicameral cycle O is not a “best” cycle, it suffices only 
(type 1 bicameral cycle for example). There may 
be other far better cycles for cycle cancellation, but it takes 
too long to compute. Secondly, a type-1 bicameral cycle Oi 
must satisfy an additional constraint 0 < c(Oi) < Copt, 
and a type -1 bicameral cycle O2 must satisfy —Copt < 
c(0) < 0. This makes the definition of bicameral cycle more 
complicated, but it is essential. Since without the constraint, 
the ratio of the algorithm will become (l + a, l-f^). In this 
case, the cost of the solution resulting from the algorithm 
could be very large when a is a small number, say a = -A 
(see Figure [ 1 ] for example). 

Because of the complicated constraints on type-1 and type- 
2 bicameral cycles, and the fact that the residual graph de¬ 
fined in this paper allows both negative cost and delay on 
edges, previous techniques as in lanani] are not suitable 
for computing a bicameral cycle. Anyhow, this paper figures 
out a constructive method for computing bicameral cycles, 
which will be shown in Section 4. 

3.2 Proof of Approximation Ratio and Time 
Complexity 

Assume that the algorithm terminates in / iterations. 
Without loss of generality, assume that / > 1. Let Oi be 
a bicameral cycle in the ith iteration, and Di and Ci be 
the delay and cost of the current solution, respectively. Let 

ADi = D — Di, ACi = Copt — Ci, and Vi — The 

1“^ 

ratio of Algorithm [1] is as below: 

Lemma 11. If the kRSP problem is feasible, Algorithm\^ 
outputs a solution with delay strictly bounded by D and cost 


bounded by 2 Copt- 

Proof. Algorithm[T]terminates only when the delay con¬ 
straint is satisfied. Then from Lemma the lemma is obvi¬ 
ously true for the delay constraint part. For the cost part, 
according to the dehnition of bicameral cycles, the cost aug¬ 
mentation in the last iteration, i.e., the /th iteration, is at 
most Copt- Then we need only to show that C/_i < Copt- 
Below is the detailed proof by using mathematical induction. 
According to Lemma [ 5 ] for the cost Co and delay Dq of the 
first phase, there exist 0 < a < 2, such that Co = aCoPT 
and Dq = (2 — a)D- Then since the algorithm didn’t termi¬ 
nate in the hrst iteration, we have Dq > D and Co < Copt- 
By induction, Ci < Copt holds. It remains only to prove 
Ci+i < Copt for the case Ci+i > Ci. That is, the com¬ 
puted bicameral cycle Oi suffices c(Oi) > 0. According to 
the definition of type -1 bicameral cycles, we have 

d{Oi) _ Dj+i — Dj ^ ADj _ D — Dj , . 

c{Oi) ^ Ci+i -Ci - 'Aa ^ Copt - Ci ' 

Before the /th iteration the delay of the solution is always 
larger than D, so Di+i > D and Di > D both hold. Thus, 
we have: 

Q> Di+i_ - Di> D - Di- (2) 

Combining Inequality m and Inequality © yields Ci+i — 
Ci < Copt — Ci, so Ci+i < Copt- This completes the 
proof. □ 

The statement of this lemma is only for the case that the 
IcRSP problem is feasible. When the instance for fcRSP is 
infeasible, i.e., there does not exist k disjoint paths satis¬ 
fying the delay constraint D, Algorithm [T] will detect the 
infeasibility and return “infeasible” in Step 2(a). 

Before presenting the time complexity analysis of Algo¬ 
rithm [U we would like first to investigate some properties of 
bicameral cycles, as in the following lemma: 

Lemma 12. For the ith iteration, i < f, at least one of 
the following two cases holds: 

1. ri+i = n and ADi+i < ADi ; 

2. ri+i > n- 

Proof. Recall that Oi is the bicameral cycle computed 
in the ith iteration. If Oi is a type-0 bicameral cycle, then 
d{Oi) < 0 and —Copt < c{Oi) < 0, or d{Oi) < 0 and 
—Copt < c(Oi) < 0 hold. For these two cases of type-0 
bicameral cycle. Clause 2 obviously holds. So we need only 
to consider type -1 and type -2 cycles, which are as below: 

1. d{Oi) < 0 and 0 < ciOi) < Copt 

According Definition 1101 < Ti holds. Then we 

have 

d{Oi) < ri - c{Oi). (3) 

After the cycle cancellation wrt Oi, ri+i = Ac‘-c(o’j ' 
Combining this inequality with Inequality © yields 


^ ADi - n - c(Oi) _ ACi - c(Oi) 
- ACi - c(Oi) ■ ACi - c(Oi) 


= ri. (4) 


Since d{Oi) < 0, ADi+\ < ADi holds. So if the two 
sides of Inequality © are equal. Clause 1 holds; oth¬ 
erwise Clause 2 holds. 













(0, (D + 1) >1= Cop^ — e) 



_Basic output . Optimal solution _Algorithm’s option 

Figure 1: An example for execution of Algorithm [T1 without consider the constraint on the cost: sabct and st 
is the solution resulted from the simple approximation algorithm in |17L [9] : sat and st would the output of 
the algorithm, with cost Copt * {D + 1) — e and delay 0; while the optimal solution is sabt and st, with cost 
Copt and delay D. 


2. —Copt < c{Oi) < 0 and d{Oi) > 0 

The proof of the second case is similar to the first case. 
Oi is with maximum according to Definition (TO] 
so > Ti holds, and hence 

d{Oi) < n ■ c{Oi). (5) 

Then combining the above inequality with Vi+i = ‘^^'Zt(o') 
we have 

^ A.Di — n ■ c{Oi) ^ ACi - c{Oi) 

ACi - c(Oi) ^ ACi - ciOi) ~ 

Therefore, Clause 2 always holds for case 2. This com¬ 
pletes the proof. □ 

Let tbc be the time of computing a bicameral cycle. We 
now consider the time complexity of the algorithm. The key 
observation is that Vi is an arbitrary number. In fact, there 
are at most ^ c(e) * d{e) different values for r. That is 
because for the formula ri = the value of ADi must 

be an integer between 0 and — a-nd so is ACi. Then 

according to Lemma 1121 the algorithm computes at most 
^ c(e) > 1 = d{e) bicameral cycles to decrease the value of 
Ti. On the other hand, the algorithm computes at most 
1^1 * X] ‘-(fi) * X] d{e) bicameral cycles to decrease the delay 
sum of the k disjoint paths (Note that ADi > ADi+i may 
hold). Therefore, we have the following lemma: 

Lemma 13. The time complexity of our alaorithm is 0(\D\* 
Y,c[e)*Y,d{e)*tbZ). 

The time complexity of our algorithm can be analyzed bet¬ 
ter. However, we omit the detailed and sophisticated anal¬ 
ysis, because the length limit, and the focus of this paper 
should be on the approximation ratio. 

4. COMPUTING BICAMERAL CYCLE Oi 

This section will show how to compute a bicameral cycle 
Oi. This task is not easy, because almost every task involv¬ 
ing bicriteria negative cycle is AP-hard. The key idea of our 


algorithm is to construct two auxiliary graphs (B) and 
such that every cycle containing v with total cost 
between 0 and P in G is corresponding to a cycle in HfZ{B) 
while every cycle containing v with cost between —B and 0 
corresponds to a cycle in HT{B), where P is a given bound 
on cost and u is a vertex of G. Then by computing the cycles 
in Hi^fB) and H~{B) for every v £ G and every necessary 
P, we could find in G bicameral cycles if there exists any. 
The construction of auxiliary graph will be given in Subsec¬ 
tion 4.1 and the other parts of the algorithm shall be given 
in Subsection 4.2. 

4.1 Construction of Auxiliary Graph p„(P) 

The algorithm of constructing the auxiliary graph (P) 
and HT{B) is inspired by the method of modeling a shallow- 
light spanning tree path subject to multiple constraints [S]. 
The full layout of the construction (P) is as shown in Al¬ 
gorithmic] (An example of such a construction is as depicted 
in Figure [2|). 

4.2 Computation of Bicameral Cycle Oi 

This subsection shall show how to use LP-rounding method 
to compute cycles in HfZ(B) and HT{B), and obtain a bi¬ 
cameral cycle in G. We start from the following linear pro¬ 
gramming formula: 

min ^ c(e)a:(e) ( 6 ) 

subject to 

E -- E x,= 0 'iv £V{H+{B)orHT{B)) 
d{e)x{e) < AD 

Unlike other linear programming formula for shortest paths, 
spanning trees, or Steiner trees, 0 < x{e) < 1 is not neces¬ 
sary for our formula, because our goal is only to compute a 
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Figure 2: Construction of auxiliary graph H with cost constraint B — 6. (a) graph G; (b) graph G wrt path 
sxyzt-, (c) auxiliary graph H. 












Algorithm 2 Construction of auxiliary graph 

Input: Residual graph G with edge cost c : e —>■ Zq" and 

edge delay d : e —>■ Zj, a specified vertex v £ V, a, given 

cost constraint B £ Z"*"; 

Output: Auxiliary graph 

1. For every vertex vi of V do 

Add B + 1 vertices vf,..., vf to (B) ; 

2. For every edge e = {vj,vi) £ E do 

(a) if c(e) > 0, add to (B) the edges 

(vj, , ■ ■ ■, each of which is 

with delay d(e); 

(b) if c(e) < 0, add to Bl^(B) the edges 

(vf, , ■ ■ ■, each of which 

is with delay d(e). 

3. For each vertex u' of H^(B) that corresponding to 
the specified vertex v £ G do 

Add edge e(n*, u°) to (B) with delay zero. /* To 
construct H~{B) is to add edge e(u', v^). */ 


bicameral cycle. Below we shall discuss Bl^ (B) only, since 
the case for H^{B) is similar. 

Lemma 14. A solution to LP m exactly corresponds to 
a set of cycles with cost between B and —B in residual graph 
G. 

A solution to LP (O is apparently corresponding to a set of 
fractional cycles of Hif[B) (or Hp{B)). So we need only to 
show that a cycle in (B) corresponds to a set of cycles in 
G. In fact, according to the construction of {B) , we have 
the following lemma, from which the correctness of Lemma 
HU can be immediately obtained. 

Lemma 15. A cycle in (B) corresponds to a set of 
eyeles in G, each of which is with cost between B and —B. 
Conversely, a cycle in G, containing v and with cost between 
0 and B. corresponds to a cycle in Hff{B). 

Proof. Assume that O is a cycle in Hff{B). According 
to the rules of adding the edges of IPf{B), every edge in 
Hff [B) either corresponds to an edge in G or is between 
duplications of an identical vertex of G. So O corresponds 
to a closed walk in G, and hence it corresponds to a set of 
cycles of G. The cost bound of the set of cycles follows from 
the construction of Hff{B) similarly. 

Conversely, let O' be a cycle in G, with cost between 0 
and B, i.e., 0 < c{0') < B. Start from € Hif(B), a 
vertex corresponding to u £ O', according to the edges on 
O', we can find a path from to \ since c(O') < B. 
The path plus the edge exactly compose a cycle 

in H+{B). □ 

Theorem 16. Let cycles 0\,...,0z be the cycles in G 
corresponding to x{v,B*) for every v £ G, where B* = 
c(Pi*) and Pf,..., Pf is an optimal solution to the 
kRSP problem. Then there must be a bicameral cycle in 
Oi,..., Oz if the kRSP problem is feasible. 


Proof. If there exist c(Oi) < 0 and d{Oi) < 0, or d{Oi) < 

0 and c{Oi) < 0, then Oi is type-0 bicameral cycle. So with¬ 
out loss of generality, we can assume that either c(Oi} > 0 
or d{Oi) > 0 holds. Let B*) be an optimal solution 

to LP (|51) against LTf{B). 

Let Pi,... ,Pk be the current solution against fcRSP. Let 
{Oi*,..., OI} = {Pi ,. .., Pfc } © (PT, .. .,T\}. Then appar¬ 
ently 

- E ) ^ - E c(P) < c(0*) < ^ c( ). 

i=l z = l i=l 

That is, cycle O* with d(0*) < 0 and c(0*) > 0 can find 
its corresponding cycle in U„ggPj'(P*). Without loss of 

generality assume that O* with minimum among all 

other cycles of {Oi,..., O^} with negative delay, and its cor¬ 
responding cycle is in (B*). Then we have || < 

I . So in the set of cycles of G that corresponding to 

{v, B*), there must exist a cycle Oi with < t[o*] 

and d{Oi) < 0, or j and c(02) < 0. The case 

for O* with d{Oi) > 0, c(0*) < 0 and maximum is 

similar in U^ggP,7(P*), and hence we can find O'l and O 2 
accordingly. Then from the definition of bicameral cycle, at 
least one cycle of {Oi, Oi, O 2 , O 2 } is a bicameral cycle. This 
completes the proof. □ 

Following Theorem ll6l the main steps of our algorithm roughly 
proceeds as: 

1. Construct (B) and Hf{B) for all n £ G and every 
integer 0 < B < Eeeo ^(e); 

2. Solve LP (0 against every Hf'{B) and Hp{B), and 
collect optimal solutions; 

3. Compute a bicameral cycle among all the cycles re¬ 
leased from all the solutions. 

In general, the above algorithm follows to the LP-rounding 
algorithm framework. Let’s revise the algorithm as well as 
its proofs following the traditional line of analyzing a LP- 
rounding based approximation algorithm. Clearly, the algo¬ 
rithm solves LP formulas, and rounds Xe to 1 for every edge 
e that belongs to a computed bicameral cycle. Then the core 
task in the analysis is to show that the cycles corresponding 
to the solution of LP always contains a bicameral cycle, 
if the fcRSP problem is feasible. This task has actual been 
done in Theorem m 

The detailed algorithm is as in Algorithm [S] 

Theorem 17. Algorithm 0 correctly computes a bicam¬ 
eral cycle in time the = 0{n‘^'^GopT^c{e)L). 

Proof. Following Algorithm |3l the resulting cycle satis¬ 
fies the definition of bicameral cycle. So it remains only to 
show the time complexity of the algorithm. Step 1 of the al¬ 
gorithm constructs 0{V{Ht {B))*Y} c(e)) = 0{nGoPT c(e)) 
auxiliary graphs, each takes 0{V{H^ (b))^'^L) = 0{n^'^0^%rpL) 
to solve the linear programming formula, where L is the 
maximum length of the input [15]. Other steps of the algo¬ 
rithm take trivial time comparing to Step 1. Therefore, the 
time complexity of the algorithm is 0{n'^'^G o%t 'FZ c(e)Z/). □ 













Algorithm 3 Computation of bicameral cycle O. 

Input: Graph G; 

Output: A bicameral cycle O in G. 

1. For = 1 to ^ c(e) do 

(a) For each v £ G do 

i. Construct (B); 

ii. Solve LP ([6| against (B), and obtain two 
optimal solution x'^{'v,B); 

iii. Release the set of cycles 0'’"(u, B) in G, which 
corresponds to x^i'>^yB)- 

iv. If there exists O £ 0~^(v,B) with c(0) < 0 
and d(0) < 0 or c(0) < 0 and d(0) < 0, 
then terminate; 

/* O is a type-0 bicameral cycle. */ 

(b) For each v £ G do 

i. Construct H~{B)-, 

ii. Obtain the set of cycles 0~(v, B) , and check 
if it contains type-0 bicameral cycles; 

/* O is a type-0 bicameral cycle. */ 

2. Compute Oi with minimum and d(Oi) < 0, and 

c(_Cyi J 

O 2 with minimum and c( 02 ) < 0 among all the 

cycles of U U (0''"(u, B) U 0~{v, B)); 

V B 

3- If l^l<if§f}|then 

return Oi; /* Oi is a type-1 bicameral cycle. */ 
Else return O 2 . /* O 2 is a type-2 bicameral cycle. */ 


The above time complexity is terrible, and can be signif¬ 
icantly improved by some sophisticated algorithm design 
techniques. For example, construction of auxiliary graphs 
for all B = 1 to 5]] c(e) is not necessary. Binary search can 
be applied here to find B*, and reduce the number of the 
auxiliary graphs constructed. Due to the length limit, we 
omit the details here. After all, the core of this paper is to 
show that the fcRSP problem admits an approximation ratio 
of (1 + e, 2 -I- e) for any e > 0. 

5. CONCLUSION 

For any constant e > 0, this paper first gave a polyno¬ 
mial time approximation algorithm with bifactor approxi¬ 
mation ratio (1 -I- e, 2 -|- e), based on improving an approxi¬ 
mate solution with bifactor ratio ( 2 , 2 ) via using bicameral 
cycles in the cycle cancellation method. Then this paper 
presented a constructive method for computing a bicameral 
cycle, by constructing an auxiliary graph and innovatively 
employing LP-rounding technique therein. To the best of 
our knowledge, our algorithm is the first constant factor 
approximation algorithm that computes a solution almost 
strictly obeying the delay constraint. We are now investi¬ 
gating the inapproximability of the fcRSP problem. 
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